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I. INTRODUCTION 

The century-old problem of describing a gravitation- 
ally bound system in an expanding universe in the frame- 
set of General Relativity has seen many attempts to find 
a solution. Despite its apparent simplicity, a full under- 
standing of the mechanisms involved when general and 
realistic systems are considered has yet to be found. This 
can be seen, for instance, by taking a look at the vast 
and often contradictory literature on one of the oldest 
proposed models, and possibly what could be considered 
the second simplest scenario after Schwarzschild-de Sit- 
ter, namely the McVittie solution [ ]. 

Our current understanding of the history of the uni- 
verse, built over an increasing amount of data from ac- 
curate measurements such as the CMB and supernovae 
distances, tells us that we are most likely living in a ho- 
mogeneous and isotropic universe which is undergoing 
a phase of accelerated expansion. The ACDM model, 
our best fit to observations so far, is far from being a 
complete model, given that it does not properly include 
baryons and that we lack a reliable understanding of star 
formation and feedback [2]. Therefore, it is evident that 
studying how bound systems feel the expansion of the 
universe, starting from collapsing star-forming matter up 
to galaxy superclusters at recent times, is of vital impor- 
tance for a better understanding of the cosmos. 

The first step in studying bound systems interacting 
with an expanding background is to consider a related 
problem: the formation and evolution of black holes in an 
expanding universe. Some of the simplest metrics which 
display these features belong to the Kustaanheimo-Qvist 
class of solutions to Einstein's equations for a comoving 
shear-free perfect fluid in a spherically symmetric config- 
uration [ ]. The oldest and perhaps most famous mem- 
ber of this class is the McVittie solution [ ] . Throughout 
the years it has been studied in this context, either in 
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its original form (see Refs. [G, 7] and references therein) 
or in a generalized version with a time-dependent mass 
[8-11] 

In this work we consider the original McVittie solution 
to Einstein's equations, which together with its physical 
interpretation has been debated for almost eighty years. 
Recently, there has been considerable advance towards 
understanding such a metric; for instance, after a long 
debate it has been proved [ 1 2] that the central object 
satisfies necessary and sufficient conditions to be charac- 
terized as a black hole, provided that the line element 
asymptotes to Schwarzschild-de Sitter at temporal infin- 
ity. In the same work, it has also been established that 
the metric possesses a singularity which lies in the past 
of every causal trajectory, the McVittie big bang. More- 
over these results have been shown to be valid for a gen- 
eralized version of the McVittie metric at least for some 
time-dependent masses [II]. 

The causal structure of the McVittie spacetime is one 
of the crucial points that have been debated in the lit- 
erature. Its analysis led the authors of Ref. [12] to con- 
clude that the internal apparent horizon present in the 
solution does asymptote to a black hole horizon. At the 
same time, following a similar path and using detailed 
numerical integrations of the light curves, the authors of 
Ref. [l-'S] found that the inner horizon ends up separating 
the boundary of the spacetime into two sections, a black 
hole horizon in its future and a white hole one in the 
past. In this brief work we want to focus on this issue 
to clarify a peculiar aspect of the McVittie solution. We 
will show in fact that the presence or absence of the white 
hole part, as discussed in Ref. ]1.H], crucially depends on 
the choice of the function describing the expansion. Inci- 
dentally, this means that the causal structures presented 
in the works cited above are not in contrast, rather they 
appear due to the particular choices made in each work. 

The aim of the theorem presented in the following sec- 
tions is to allow one to recognize the asymptotics, and 
thus the possible embeddings of the spacetime once the 
expansion function is chosen. Although we succeeded in 
describing a wide class of possible solutions, the theorem 
falls short of giving a prediction for a specific class of 
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expansion functions that decay exponentially in a very 
peculiar way (as will be explained in detail in section 
IV A) for which a case by case study is necessary. On the 
bright side, all realistic models of expansion nicely fit into 
the two groups of functions for which the theorem is able 
to identify the corresponding spacetime causal structure. 

After a brief summary of the main characteristics of 
the McVittie solution in section II, in section III we per- 
form an analysis of the inner apparent horizon, around 
which the behavior of geodesies decides the asymptotes 
of the spacetime. In section IV we present the main re- 
sult of this paper, a theorem that allows us to find out 
whether the inner horizon is an accumulation point from 
both above and below or just on one side. Finally we con- 
clude in sections V and VI with an instructive example 
— ACDM — and some remarks. 

Throughout the paper we use the convention of repre- 
senting the derivative with respect to the time coordinate 
t with a dot, and with respect to the radial coordinate r 
with a prime. 



II. MCVITTIE CAUSAL STRUCTURE 

The McVittie metric with constant mass parameter 
771 > can be written in the form [II, 12] 



ds^ 



H\^)df~'^drdt+^+r^dn\ (I) 



where R{r) — y^l — ^'"/r, and H(t) = °-/a is the Hubble 
parameter. H is assumed to be a smooth function with 
the following properties: 



Um H{t) = Ho , 

t^OO 

>Ho>0, 



3^/3m 
H{t)<0, t>0. 



(2a) 
(2b) 
(2c) 



Drawing the Penrose diagram of a fully dynamical met- 
ric such as (I), which satisfies certain conditions as de- 
scribed in Ref. [14], reduces to integrating geodesies and 
is thus mostly a numerical effort. As already pointed out, 
we will make some remarks on previous works on this 
topic, and for an easier comparison we will then adopt 
the choices made in Ref. [13]. We thus choose the sin- 
gular surface at r = 2m, which lays in the past of every 
event, as our reference point, identifying it with a hor- 
izontal line in the conformal diagram. All ingoing and 
outgoing geodesies may be interpreted as departing from 
this surface, and thus every event in spacetime is selected 
by two coordinates. These are the time at the singular- 
ity of the ingoing and outgoing geodesies meeting at the 
event. We take the function transforming from times to 
coordinates in the causal diagram to be equation (35) in 
Ref. ] ]. 

The behavior of the outgoing geodesies is well under- 
stood, and their integration does not present difficulties. 



On the other hand, due to the fact that the apparent hori- 
zons of the McVittie metric are anti-trapping surfaces, 
and thus influence the behavior of ingoing geodesies, the 
latter are to be treated more carefully and present the 
only source of possible confusion. In particular, the pres- 
ence of an accumulation point for the ingoing geodesies 
makes numerical analysis and integration challenging. 
Therefore, in what follows we focus on such geodesies 
represented by solutions of the differential equation for 
the function r : (io,oo) — (2m, oo) ]i2] 



r{t) = R{r) [rH{t)~R{r)]=X{t,r), 



with initial condition 



riU) = ro > 2m , 



(3) 



(4) 



which ensures that the coordinates in (I) describe the 
physical space above the singularity. Note that property 
(2c) imphes X < 0. 

By integrating backward in time the ingoing geodesies, 
we can define a first time to, the time at which each 
geodesic leaves the singularity at r = 2m, taken to be the 
reference for the conformal diagram. Correspondingly we 
will refer to the geodesic starting at to as rtg{t), so that 

lim Tfg (t) = 2m . 

Let us denote f{t,r) — ^j^j^', the apparent horizon % 
is the locus of points in the (t, r) plane in which f[t, r) — 
0. In other words, the horizon is the set of points where 
X{t, r) = above the singularity. An example of horizons 
for the McVittie metric for a specific choice of scale factor 
is depicted in figure I. 



2m 



r , 




Figure 1 . An example of the McVittie apparent horizons plot- 
ted in spacetime coordinates. Some ingoing (r-in) and outgoing 
geodesies (rout) are also shown. The singularity at r = 2m is 
represented by a thick dotted Hne. 

It can be easily proved by convexity arguments that 
for each fixed value of t, / has up to two real roots for 
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positive values of r, and the roots always belong to the 
interval (2m, oo). Moreover, if f{t,r) has at least one 
real root, then f{t,r) has two real roots for every t > i. 
It can also be noted that when the two roots are distinct, 
and thus have multiplicity one, one has f'{t,r) ^ at 
the root. 

We call inner horizon "H- the locus defined by the 
smallest real root of f{t,r) for each t for which it ex- 
ists. Calling iinin the smallest value of t for which 
f{t, r) admits at least one real root, the inner horizon 
is then represented by the curve (t,r_(t)), where 
T-^ ■ [tminjOo) — > (2771, oo) Is & smooth function. We will 
call r_ the inner horizon function, or, when there is no 
room for confusion with the set H^, just inner horizon 
for short. All points of correspond to points in which 
the geodesic fiow is horizontal in the (t,r) plane, as the 
tangent to the flow is given by the value of X{t, r). 

A useful characteristic of the inner horizon function 
is to have a negative slope for all times for which it is 
defined. Explicitly, by taking the gradient of f{t, r) along 
the curve, the slope reads 



-it) 



H{t)-R'{r^{t)) 



(5) 



provided that f{t,r_{t)) = H{t) - i?'(r_(t)) ^ 0. The 
numerator of (5) is clearly positive following (2c), and 
it is easy to show that the denominator /'(t, r_(i)) < 
0. In fact, noticing that f{t,2m) = 2mH{t) > 0, by 
continuity for any time for which r_(t) is a simple root, 
it follows that f{t,r^{t) — 6) > 0, for 6 > 0, ensuring 
/' negative on the inner horizon. This also implies that 
the ingoing geodesies increase monotonically in the inner 
region 2m < r < r_. 



Picard-Lindelof theorem, which means that the ingoing 
geodesies at each event above the singularity are unique. 
Then, the image of an event by the first-time projector 
$, if it exists, is a unique first time tQ. 



III. INNER HORIZON ASYMPTOTIC 
BEHAVIOR 

Why are we interested in the image of the geodesic 
fiow projected on the singularity? The main reason for 
that can be understood by looking at two cases consid- 
ered in the literature. In Ref. [12] one can see that the 
image oi %- under the projector $ defined in section 
II A is unbounded, and the asymptotic analysis shows 
the presence of a black hole (see figure 2) . On the other 
hand, the analysis done in Ref. [i:!], where the projection 
of the inner horizon is bounded, finds that the space- 
time presents a bifurcation two-sphere that separates an 
asymptotic black hole and a white hole (see figure 4). We 
agree with the results presented by both groups, and we 
will show in what follows that the discrepancies are due 
to a different choice for the function describing the expan- 
sion of spacetime. Moreover, we want to show that the 
connection between the projected flow and the asymp- 
totic behavior is what allows one to distinguish between 
the two possibilities of inevitably flnding a white hole ac- 
companying the black hole or not. So the problem we 
are facing can be restated as follows: is the image of T-L- 
under <i> bounded? A negative answer would imply the 
presence of a black hole alone, while a positive one would 
change the structure at inflnity and produce a white hole 
as well. 



A. The first-time projector $ 

To prepare for our main result, we are going to intro- 
duce a function which will play a crucial role in the rest 
of this work, the first-time projector $. This function, 
or better the image of the inner horizon under it, will 
contain the information on the causal structure of the 
spacetime considered. 

Let us call M. the set of events covered by our coordi- 
nates (t,r), so that M C (0,c») x (2m, c»), and define 
the application $ : — ?> T, where I C (0, oo), which 
associates to each element of M the first time of the 
ingoing geodesic rtp(i) which passes by that point. In 
other words, the image of <i> is the image of the geodesic 
fiow of over the singularity; it is the set of values of 
to one can reach by integrating geodesies, starting at any 
event in Ai back to the singularity. 

It is easy to show that $ is well defined, in the sense 
that the image of each coordinate pair (t, r) in its do- 
main is at most one point of the interval (or empty) ; in 
fact, it is straightforward to check that the system given 
by equations (3) and (4) satisfies the hypotheses of the 




r = 2m 



Figure 2. An extension of McVittie spacetime with an un- 
bounded image of the inner horizon over the singularity. In- 
going geodesies are represented by lines inclined 45° to the 
left. The left part of the graph shows the appearance of a 
black hole whose horizon is given by r_ in the limit t ^ oo. 

To answer the question of boundedness, we will need 
some preliminary results in order to reduce the prob- 
lem to a simpler one. First, let us define the limit- 
ing values at temporal infinity for r_(t) and f{t,r) as 
Too = limf^oo r-{t) and fcoir) = limt^oo f(t,r). Then, 

2m < Too < r_ (t) , 
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and 

/oo(?'oo) = raoHo - R{roc.) = . (6) 

We note that, in the case in which the image of the 
inner horizon by the first-time projector $ is bounded, 
then there exists such that, for each r > ip, the ingo- 
ing geodesic r^- (t) that leaves the singularity at instant r 
never crosses the inner horizon r_(i) and tends to it at 
time infinity. This means that C (0,to]- 

In what follows we show some consequences of taking 
as a hypothesis the fact that "I>(?^_) is bounded. 

Proposition III.l. Let rt„(t) be an ingoing geodesic. If 
rtg satisfies rt„{t) < r^it) for all t > t^, then 

rtoit) < roo, Vt > to, (7) 

and 

lim rtg {t) = Too • (8) 

Proof. Consider the sequences t„ = to + nr, rn = 7'to(tn), 
hn = r_(i„), n G N, T > 0. By hypothesis, and by 
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Figure 3. A possible null geodesically complete extension of 
figure 2, where a time-reversed McVittie metric with a Big 
Crunch also appears [15]. 
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Figure 4. A null geodesically complete extension of a McVit- 
tie spacetime characterized by a bounded upper value for 
the first-time projector of the inner horizon. Notice 
that geodesies starting after do not reach the horizon in a fi- 
nite time. This completion shows the presence of a white hole 
and a black hole region, separated at the branching surface 6, 
as discussed in Ref. [io]. 



making use of the monotonic increase of inner geodesies, 
the two sequences satisfy the following properties: 

2to < r„ < hn , (9a) 
r-a > n> m. (9b) 

To prove (7), we consider by reductio ad ahsurdum the 
case r„ > r^, for some n. Then, as lim„_>.oo h^ = r^o, 
this would imply r„ > hm, for some m > n. But, by 
property (9a), < hm, which implies r„ > for n < 
TO, and contradicts property (9b). The case r„ = Too 
follows by noting that if it holds for some n, then r„_|_i > 
i^n > Too, by property (9b), which leads us to the previous 
case. 

Moving on to (8), note that r„ is monotonic and 
bounded, therefore it converges. Let us call lim„^oo = 
r* < roo its limit. Recalling the definition of X(t,r) in 
(3) and the mean value theorem, there exist sequences f„ 
and r„ = rtg{tn) such that 

rn+i - rn = TX{tn, r„) . (10) 

Since lim„_>.oo — ?'*, the limit of (10) is = 
TR{r*)foo{r*), then for this to hold it must be that 

r*=roo. □ 

The importance of proposition III.l is clear; in fact, 
thanks to it one only has to analyze geodesies crossing 
Too, a fixed surface, eliminating the complication of hav- 
ing to consider the time- varying inner horizon. Moreover, 
it implies that the Too surface behaves as an accumulat- 
ing point for the geodesies as we will see in the following 
corollary. 

Corollary III. 2. // rtg (t) is an ingoing geodesic, then 
for all e satisfying < e < r^o — 2to, there exists i > t^ 
such that rtg{t) = r^o — £• 

This means that, given enough time, all geodesies ei- 
ther cross Too, or reach values arbitrarily close to it. 
By Proposition III.l, every ingoing geodesic that never 
crosses the inner horizon never reaches r^o- Conversely, 
if an ingoing geodesic does traverse r^o in a finite time in- 
terval, then it eventually crosses the inner horizon. Then, 
by studying only the neighborhood of r^o , we may tell if 
geodesies do or do not cross the horizon. 

To complete our set of preliminary results, we will show 
now that geodesies that start at a later first time remain 
below (that is, at smaller values of r than) those which 
start at earlier first times. This can be stated precisely 
by the following. 

Proposition III. 3. Let rt^{t) and rt^it) be two ingoing 
geodesies. Therefore, if ti < t2, then rt^ > rt^. 

Proof. Using the fact that the ingoing geodesic passing by 
any event (t, r) of the plane is unique for t > 0, r > 2m, 
then we conclude that they cannot cross in the same re- 
gion. Moreover, as for any i5 > 0, rtg{to + S) > rtg+sito + 
S) — 2m, then by continuity rtg(t) — rtg^s{t) > for all 
t>to. □ 
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Propositions III.l and III. 3 imply that, if we prove that 
there exists one ingoing geodesic that never reaches r^o, 
then is bounded, as every geodesic which leaves 

the singularity later will never reach the inner horizon as 
well. Otherwise, is unbounded if every ingoing 

geodesic reaches r^o in a finite time. 

Using Corollary III. 2, we only need to study the 
geodesic flow in a small neighborhood below Too, as we 
know that it represents the point separating geodesies 
crossing the inner horizon or staying under it. 



Then, there exists 5^ > 0, continuous in e, such that 
lim£_>.o i5e = and 

Se > 0{z) > -(5, , 

for < z < e. We also define the curves z^{t), z^{t), 
which are solutions of 

4(0 = -X{t, Too) + [dX{t,roo) ± <Se] 4 : (H) 

with initial condition 



IV. THE MAIN RESULT 

We are finally ready to state the main result of this 
work, a theorem which allows one to find out whether 
all geodesies leaving at late times from the singularity 
are bound to cross the McVittie inner apparent horizon, 
or if there exists a time for which all geodesies leaving 
the singularity at subsequent times never reach the inner 
horizon and accumulate under it. In other words, we 
want to state here how the form of the Hubble parameter 
H{t) can select the properties of the image of the inner 
horizon under the first-time projection function $. 

Theorem IV. 1. Let be the real-valued function AH{t) — 
H{t) — Hq, and the constants A R{raa) + r^oR' {^00), 
B — R{r ao) {R' {foe) ^ Ho) ti > 0. // there exists 
(5 > such that 

F+{t,,t) = j\^^''^^e-''^".^"^'^'''AHiu)du 

diverges as t ^ 00, then the image of "H- under $ is 
unbounded. 

Analogously, if there exists S > such that 

F^it,,t) = ^*e(^+*>e--^-^'"^^(^'''^AiJ(u)du 

converges as t ^ 00, then the image of "H- under $ is 
bounded. 

Proof. Let us analyze the flow of equation (3) near r^o, 
since we have already seen in section III that geodesies 
that cross r^o are bound to reach the inner horizon in 
a finite time. Consider < e < r^o — "^rn. Let z{t) = 
Too -r{t), satisfying 

i(t) = -X(i,roo - z) , 

and 

z{ti) = ZQ, 0<ZQ<e. 

Now, the crossing of Too by an ingoing geodesic is equiv- 
alent to a change of sign in z{t). 

By the definition of the differential of X (with respect 
to the second variable), at r = Too, we may write 

i(t) = ~X{t, roo) + dX{t, r^)z + 0{z)z . 



4 (^i) ~ ^0 ' 

and verify zLit) < z{t) < z^{t) as long as |z±(i)| € 
[0,e). Therefore, if zl(t) is always positive, then z{t) 
never changes sign and the geodesic r{t) never crosses 
Too. Conversely, if z^{t) does change sign, then so does 
z{t), and consequently r{t) crosses Too. 

The next step of the proof is to solve (11), which is 
a non-homogeneous linear differential equation, formally 
solvable by the method of variation of constants. Explic- 

itly, 



dX{t,r^) = 




^R'{roo)[roo{Ho + AH{t))-Riroo)] 
+ Rir^)[Ho + AHit)-R'ir^)] , 

which, after using equation (6) and some rearrangement, 
gives us 

dX{t,r^) = [R'{ + R{r^)]AH{t) 

+ Rirao)[Ho~R'ir^)] 
= AAH{t)-B. (12) 

We note that the last constant term in (12) is 

-B = R{r^)fUroo) - r^H^ - < 0, 

as / is positive between the singularity and the inner 
horizon. The case /4)(^oo) = corresponds to the ex- 
tremal case where / has two coincident real solutions and 
the spacetime asymptotes to an extremal Schwarzschild- 
de Sitter, which we discarded by hypothesis with prop- 
erty (2b). Thus, B >0. The constant 

TTl 

A = R'{roo)roo + ^(^oo) = , ^ „2\ + ''oo^^o 

('■So^o) 

is also strictly positive. 

The inhomogeneous term X(t, Too) in (11) can be writ- 
ten as 

X{t,r^) = CAH{t), (13) 
where we define the positive constant 

C = R{r^)r^ ^ ri,Ho > . 
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Then, substituting equations (12) and (13) into (11), we 
obtain 

4 (t) = [AAH{t) -B±5,] z| (t) - CAH{t) , (14) 
whose solutions are, with the initial conditions Zj.{ti) = 



:g _ C / e(^T*^ e"^ ^^'^^'^^Ai7(u)dt 



(15) 



Next we need to find out under which conditions z^. 
change sign or not. Note that only the factor between 
brackets in (15) can be non-positive, depending on the 
values taken by the integrals 



which means that becomes negative for finite t, inde- 
pendently of the initial time to. Following the reasoning 
of changes sign, and all ingoing geodesies even- 

tually cross the inner horizon independently of the time 
at which they leave the initial singularity. □ 

Summarizing the results of the theorem, we have that 
the divergence (convergence) of the integral defining 
{F-) allows us to find curves below (above) any geodesic 
around Voo forcing them to cross (to stay below) the in- 
ner horizon. Crossing the inner horizon in a finite time 
independently on the starting point means of course that 
the geodesies that reach it leave the singularity at all 
times, while finding geodesies that never reach the in- 
ner horizon means that there exist a point in time where 
horizon-crossing geodesies accumulate. 



(BT5.)(«-t.)e-^/*:^^(^)d^Ai7(w)dt 
Ce(-^±*»)*'F^(i„t). 



In particular, the convergence of or the divergence 
of F^ in the i — > oo limit immediately tell us about the 
behavior of z^.. Given that F_{ti,t) > F^{ti,t), for all 
e > there are only three possible cases: 

(a) There exist M{ti) > and e > such that 
\imt^^F^{U,t)=M{tiy, 

(b) For all e > 0, limj^oo F^it^t) = oo; 

(c) There exist N{ti) > and e > so that 
\imt^^ Fl{ti,t) = N{U) but \imt^^Fi{ti,t) = oo 
for all e > 0. 

Case c does not respect the hypothesis of the theorem, 
rather it is the case in which the method here presented 
cannot be applied, and we will discuss it later. We can 
start then with case a in which the term between brackets 
in (15) becomes 



Ce 



M{U) 



However, as limt.^oo Ce ^*-^M{ti) — 0, there exists r > 
such that for alH > r, z^-Ce~^*M{t) > 0. This means 
that after the instant r the zl curves do not change sign 
anymore, and since they are a lower bound for z, neither 
do the curves with ti > t. It follows that the ingoing 
geodesies rtg(t), which reach Too— Zg times equal to r or 
later, do not cross the inner horizon at a finite coordinate 
time. As to < r, this gives us an upper bound to the 
image of the inner horizon under the first-time projector 
as $CH_) C (0,r]. 

In the case described by b instead, we see that there 
is no upper bound to $('H_). In fact, in this case there 
always exists a time Tf . > such that, for each ti > 0, 



A. Limits of applicability 

As we said before, case c describes expansion functions 
for which we cannot immediately apply our method to 
find the asymptotes of the spacetime. The Hubble pa- 
rameters that fall into this scenario are those that take 



the form AH{t) 



-Bt 



h{t), with B the constant defined 



in the theorem, for which the conditions 

(i) Ve > 0, Vt, > 0, the integral J^"" e-'*h{t)dt con- 
verges, 

(ii) yti > 0, the integral h{t)dt diverges, 

are both satisfied. These choices describe the "blind spot" 
of the method we presented. Although there is an infi- 
nite number of functions which may be constructed with 
these properties, the fine-tuning required by the expo- 
nential part means that such functions do not constitute 
a significant fraction of physically relevant cases, as we 
will see in an example in the next section. 



V. EXAMPLE: ACDM EXPANSION 

Let us illustrate the method described by Theorem 
IV. 1 by applying it to an example, which will also help 
showing that most physically relevant expansion func- 
tions fall within the purview of the method. We choose 
the Hubble parameter used in ACDM models where dark 
matter and dark energy are the main components of the 
energy budget, a good approximation to the universe as 
we see it today, and also for the description of large-scale 
structure formation. In this case, H can be described by 
[13, 16] 



H(t) = iJoCoth ( ^Hot 



(16) 



Keeping in mind that we want to analyze the asymp- 
totes of the spacetime, we may expand (16) at late times 
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and rewrite it as 

H{t) = Ho + 2Hoe-^"°* + 0{e'^"°') , 

which corresponds to AH = 2iJoe~^^°* [l + 0(e"^^°*)] . 
To prove that in this case the image of H- is bounded, 
we need F„ (i,; , t) to converge as per case a. On the other 
hand, if F^{ti, t) diverges, then we will have proved that 
the image of is unbounded by falling into case b of 
the theorem. 

With this form for AH, the function F_ is given by 
F^{t,,t) = 

Since the last factor is bounded between the two positive 
values 

e ^ < e -"i < 1, 

the convergence of is determined only by the integral 



therefore, F_ will converge if there exists (5 > such that 
B - 3Ho + 6 <0. This is true if and only if 

B-3Ho<0. 

Following the same reasoning, we find that F+ diverges 
if and only if 

B-3Ho>0. 

Inserting the definition for B, we notice that the rele- 
vant parameter for our discussion is given by 



V ■ 



B 
3i^o 



1 



R{ro 



Ho 



1 



1 



(17) 



Therefore, if 77 < then F^ converges, and we have 
a limiting first time. If 77 > then F^ diverges, and all 
ingoing geodesies cross r_. If 77 — 0, then F+ converges 
and F^ diverges, satisfying both properties i and ii in 
section IV A, and therefore leaving us with no knowledge 
about the ultimate fate of geodesies in the spacetime. 

To calculate the values 77 may assume we use the 
fact that the constant B can be explicitly calculated in 
terms of r^o ■ For McVittie metrics which asymptote to 
Schwarzschild-de Sitter, the value of r^o is given by [17] 



cos 



1 



(^SV^mHo 



(18) 



It is easy to show that, once we insert (18) in (17), all 
dependence with respect to the free parameters 771 and 
Hq is expressed in terms of the product A = mHo, i.e., i] 



M>0 \ 

diverges (b) 


/ n = (c) 




F_ converges (a) 



0.0755 

k = tuHq 



Figure 5. Convergence criterion applied to ACDM (16). Val- 
ues of A = mHo situated to the left of the root correspond to 
inner horizons which are always reached in finite time. Values 
situated to the right correspond to metrics with horizons that 
are only reachable by geodesies which start before a critical 
first time. 



is constant along hyperbolas in the m-Ho plane. More- 
over, by noting that 771 > and by making use of (2b), 
we find that the region of the parameter space available 
corresponds to < A < We plot the values taken 

by rj within this region in figure 5. 

In other words, in the region of the parameter space 
where 77 < the causal diagram will be like the one in 
figure 4, whereas where 77 > the causal diagram will 
resemble the one from figure 2. In this model, 77 has a 
root for A = 2 — ~ 0.0755. Only at this very curve 
our criterion fails to determine whether the image of 
is bounded, therefore in this case nothing can be said as 
to the final fate of geodesies. 



VI. CONCLUSIONS 

In this paper we developed a working method to deter- 
mine the causal structure of McVittie metrics using the 
information contained in the Hubble parameter. One of 
two outcomes is possible: either geodesies leaving the 
singularity at all late times cross the inner horizon, cor- 
responding to a case in which <i>(H_) is unbounded, or 
there exists a time for which any geodesic leaving the sin- 
gularity after this upper bound never reaches the inner 
horizon, corresponding to a bounded <I'('H-). The causal 
structure that results from the first case is shown in fig- 
ures 2 and 3, where the inner horizon at time infinity 
identifies with a black hole event horizon. The second 
case results in a causal structure as in figure 4, where a 
bifurcating surface appears, splitting the boundary into 
a black hole in its future and a white hole in its past. 

We showed that we can reduce the problem to the anal- 
ysis of whether ingoing geodesies cross a fixed threshold 
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given by the position of the inner horizon at future infin- 
ity. By using the Hnearized form of the geodesic equation 
we constructed majorant and minorant functions which 
can be treated analytically, and whose behavior near the 
threshold, due to the fact that they shepherd the geodesic 
between them into the same region, can be used to trace 
the solution of the full geodesic equation. The method 
is inconclusive only if the majorant crosses the threshold 
and the minorant does not, a situation which does not 
fix the behavior of the geodesic in between. 

We applied this formalism, which is stated formally 
in Theorem IV. 1, to analyze the causal structure of the 
McVittie metric when the expansion factor is given by a 
ACDM model. The simplicity of this example allowed us 
to cover the full spectrum of accessible values for the two 
independent parameters left in the metric, namely the 
black hole mass m and the constant Hq, the asymptotic 
value of the Hubble parameter at time infinity. We found 
that both causal structures are possible, depending on 
the values of these parameters, and that the method we 
developed only fails at the curve given by 77 = in the 
allowed portion of the two-dimensional parameter space 
for m and Hq. 

The aim of this work is not to discuss the meaning of 



the asymptotics of McVittie spacetime, especially since 
this has been done by various other authors (see for in- 
stance Refs. [12, 13, 15]). Here the focus has been on 
finding a way to resolve the confusion present in the lit- 
erature, and to be able to tell apart completely different 
physical setups that are generated from the same metric 
via different choices for the Hubble parameter. To clarify 
this important point, we used examples borrowed from 
the literature, in particular identifying the two possible 
scenarios with the two cases studied in Refs. [12] and 
[1.!] as shown in figures 2 and 4, even though the results 
discussed here are fully general. With theorem IV. 1 we 
proved that the form of the expansion function is the sole 
responsible factor for the structure of the boundaries in 
McVittie spacetime. 
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